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ABSTRACT 

We perform a cross-correlation analysis of microwave data from Wilkinson Microwave Anisotropy Probe 
and photometric quasars from the Sloan Digital Sky Survey, testing for Sunyaev-Zeldovich (SZ) effect from 
quasars. A statistically significant (2.5 o) temperature decrement exists in the 41 GHz microwave band. A 
two-component fit to the cross-correlation spectrum incorporating both dust emission and SZ yields a best-fit 
y parameter of (7.0 ±3.4) x 10^^. A similar cross-correlation analysis with the Luminous Red Galaxy sample 
from Sloan gives a best-fit y parameter of (5.3 ±2.5) x 10^^. We discuss the possible physical origin of these 
signals, which is likely a combination of SZ effects from quasars and from galaxy clusters. Both the Planck 
Surveyor satellite and current ground-based arcminute-resolution microwave experiments will detect this signal 
with a higher statistical significance. 

Subject headings: Cosmologyxosmic microwave background-galaxies:active-submillimeter: general 



1. INTRODUCTION 

Secondary anisotropies in the cosmic microwave back- 
ground (CMB) provide a number of opportunities for prob- 
ing the growth of structure in the Universe. Numerous recent 
and current experiments aim to detect these effects. The ther- 
mal Sunyaev-Zeldovich effect (TSZ) (Sunyaev & Zeldovich 
1972), which is the inverse Compton scattering of cosmic 
microwave background (CMB) photons due to hot electrons, 
creates a spectral distortion in the CMB and is the dominant 
secondary effect contributing to arcminute-scale anisotropies. 
The TSZ effect arises from the random thermal motion of the 
electrons. If there is a finite velocity of the cluster in the 
CMB frame there will be an additional Doppler term. This 
Doppler anisotropy is called the kinetic SZ (KSZ) distortion. 
The TSZ effect serves as a powerful tool for probing accumu- 
lations of hot gas (see Carlstrom, Holder, & Reese 2002 for 
a review). Most of the TSZ signal comes from galaxy clus- 
ters, which contain the majority of the thermal energy of the 
Universe. State-of-the-art SZ surveys like the Atacama Cos- 
mology Telescope (ACT') (Hincks et al. 2009) and the South 
Pole Telescope (SPT^) (Staniszewski et al. 2009) aim to de- 
tect large numbers of clusters in blind surveys via their SZ 
signatures. 

A number of other astrophysical processes will also create 
both TSZ and KSZ distortions. These include KSZ distortions 
from peculiar velocities during reionization (McQuinn et al. 
2005, Illiev et al. 2006), both thermal and kinetic distortions 
from supernova-driven galactic winds (White, Hernquist, & 
Springel 2002; Majumdar, Nath, & Chiba 2001), thermal ef- 
fect from black hole-seeded proto-galaxies (Aghanim, Bal- 
land, & Silk 2000; de Zotti et al. 2004, Rosa-Gonz'alez et al. 
2004, Massardi et al. 2008), KSZ from Lyman-break galaxy 
outflows (Babich & Loeb 2007), and thermal SZ effect from 
supernovae from the first generation of stars (Oh, Cooray, 
& Kamionkowski 2003). Here we consider another generic 

' http://www.physics.princeton.edu/act/ 
^ http://pole.uchicago.edu 



signal: the TSZ distortion from hot gas surrounding active 
galactic nuclei (AGN) (e.g., Natarajan & Sigurdsson 1999; 
Yamada, Sugiyama, & Silk 1999; Lapi, Cavaliere, & De 
Zotti 2003; Platania et al. 2002; Roychowdhury, Ruszkowski, 
& Nath 2005; Chatterjee & Kosowsky 2007; Scannapieco, 
Thacker, & Couchman 2008; Chatterjee et al. 2008). The SZ 
effect from these sources provides a new observational tool to 
study the role of baryonic physics in structure formation. 

X-ray evidence of the effect of AGN feedback on struc- 
ture formation at different scales has been shown by several 
groups (e.g., Mcnamara et al. 2005; Voit & Donahue 2005; 
Sanderson, Ponman, & O' Sullivan 2006; Weinmann et al. 
2006; Schawinski et al. 2007). At the same time, study of ra- 
dio jets has illuminated the morphology and outburst history 
of AGNs (e.g., Birzan et al. 2004; Dunn, Fabian, & Taylor 
2005; Dunn & Fabian 2006). Although effects of AGN feed- 
back have been successfully detected on different scales, the 
challenge lies in quantifying the amount of feedback energy 
and its particular effects on structure formation. Adding SZ 
(thermal from now on) measurements to X-ray data will pro- 
vide stronger constraints on feedback effects than X-ray data 
alone (e.g., Scannapieco, Thacker, & Couchman 2008; Chat- 
terjee et al. 2008; Moodley et al. 2008). But the AGN signal is 
tiny, and direct detection through pointed observations with, 
for example, the Atacama Large Millimeter Array (ALMA^) 
(Lapi, Cavaliere, & De Zotti 2003; Chatterjee & Kosowsky 

2007) will be challenging. In this work we instead cross- 
correlate CMB maps with optical quasar catalogs (Chatter- 
jee & Kosowsky 2007, Scannapieco, Thacker, & Couchman 

2008) to enhance the detectability of the SZ signal. 

The cross-correlation function and its Fourier transform, 
the cross-power spectrum (Peebles 1980), have been powerful 
techniques in cosmology to probe small signals. The cross- 
correlation of CMB data sets with galaxy surveys has proven 
useful for measuring secondary anisotropies in the CMB (e.g., 
Peiris & Spergel 2000; Refregier, Spergel, & Herbig 2000; 
Fosalba, Gaztanaga, & Castander 2003; Afshordi, Loh, & 

^ http://www.alma.nrao.edu/ 
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Strauss 2003; Padmanabhan et al. 2005; Ho et al. 2008; Gi- 
annantonio et al. 2008; Hirata et al. 2004; Smith, Zahn, & 
Dore 2007; Diego, Silk, & Sliwa 2003; Cheng, Wu, & Cooray 
2006; Croft, Banday, & Hemquist 2006; Myers et al. 2004; 
Kashlinsky et al. 2009; Diego & Partridge 2010; Hernandez- 
Monteagudo et al. 2006; Ho, Dedeo, & Spergel 2009). Here 
we investigate the SZ distortion from AGN feedback by cross- 
correlating the WMAP CMB maps with a photometric quasar 
catalog selected from the Sloan Digital Sky Survey (SDSS). 

Our paper is organized as follows. In §2 we give a brief 
description of the data sets used in this work. Section 3 dis- 
cusses our methodology for estimating the amplitude of the 
SZ-quasar cross-correlation. In §4 we describe potential sys- 
tematic effects and the masks used to eliminate data which is 
likely to be dominated by systematic errors. Our results are 
discussed in §5. Finally, in §6 we give a brief summary of 
our results and discuss the prospects of detecting quasar-SZ 
cross-correlations with future experiments. 

2. DATA SETS 

The analysis presented here uses four data sets: the Wilkin- 
son Microwave Anisotropy Probe (WMAP) CMB tempera- 
ture maps, the SDSS quasar and luminous red galaxy (LRG) 
catalogs, and the NRAO VLA Sky Survey (NVSS) radio cat- 
alog. Table 1 summarizes the sky areas, objects, and redshift 
ranges for these surveys. Figure 1 displays the sky regions 
covered by the SDSS and NVSS catalogs, as well as the masks 
used in the analysis. 

2.1. WMAP Temperature Maps 

The WMAP satellite has mapped the microwave sky in 5 
frequency bands (Bennett et al. 2003) with an angular resolu- 
tion ranging from 0.88 degrees to 0.22 degrees (corresponding 
to multipole moments ranging from / = 250 to / = 900. We 
have used the WMAP 5-year foreground-reduced CMB tem- 
perature maps for Q (41 GHz), V (61 GHz), and W (94 GHz) 
bands (Gold et al. 2009). The temperature 7] assigned to pixel 
; is the temperature difference from the mean map tempera- 
ture, after dipole and foreground removal. These maps use 
the HEALPix scheme developed by Gorski et al. (2000). We 
have used the HEALpix resolution 9 maps for our analysis. 
This corresponds to a total of 3145728 pixels with each pixel 
being 47.2 square arc minutes in area. For each frequency 
band we use a weighted average of the maps from the individ- 
ual differencing assemblies (Bennett et al. 2003). 

2.2. SDSS Quasar Catalog 

SDSS has performed 5 band (ugriz) photometry (Fukugita 
et al. 1996) over roughly 10000 square degrees of sky area. 
We use the quasar catalog developed by Ho et al. (2008). The 
quasars are selected photometrically. A quasar candidate cata- 
log of UV-excess objects {{u — g) < 1 .0) with observed g mag- 
nitudes fainter than 14.5, extinction corrected g magnitudes 
brighter than 21.0, and u-band error less than 0.5 mag are 
generated (Ho et al. 2008). A catalog of DR3-UV-excess ob- 
jects is also constructed. The candidate quasar catalog is then 
matched in color space with the DR3-QSO catalog (Richards 
et al. 2006) and the DR3-UV-excess catalog to generate the 
current QSO catalog. Three additional maps were used to 
clean the data set from stellar contamination and mask out re- 
gions affected by poor seeing: (i) a point spread function map; 
(ii) a stellar density map (18.0<r< 18.5 stars smoothed with 
a 2-degree FWHM); and (iii) a stellar surface density map us- 
ing red stars with {g — r) greater than 1.4. These maps are 



TABLE 1 
Data Set SPECincATioNS 



Catalog 


Redshift 


Area 


Objects 






(deg2) 




QSO catalog 


0.08-2.82 


6039 


586435 


LRG catalog 


0.4-0.6 


6641 


911686 


NVSS catalog 


0-3.0 


27361 


1104983 



described in detail in Ho et al. (2008). The data set here uses 
SDSS PSF magnitudes for all quasars. 

2.3. SDSS Luminous Red Galaxy Catalog 

Luminous red galaxies are useful mass tracers in the Uni- 
verse. The LRGs tend to have old stellar population with 
a uniform spectral energy distribution. The characteristic 
4000 A break makes LRGs ideal objects for photo-z algo- 
rithms with redshift accuracies « 0.03 (Padmanabhan et 
al. 2005). The selection criterion, color and magnitude cuts 
for the LRG sample employed here are described in Pad- 
manabhan et al. (2005). The LRGs are chosen to be galax- 
ies that have colors consistent with an old stellar popula- 
tion, and absolute magnitudes above a certain threshold (Pad- 
manabhan et al. 2005). The magnitude cuts are given as 
'-petro < (13.6 + 2.3(^-r)+4(r-/- 0.18)), rpetro < 19.7, 
/ < (18.3 + 2{r - i) - 0.25{g - r)), and / < 20 where rpetro 
is the SDSS r-band Petrosian magnitude (Padmanabhan et al. 
2005). The LRG catalog used for the current work is drawn 
from the SDSS sample prepared by Ho et al. (2008). The 
magnitudes we use for the LRGs are model magnitudes (as 
opposed to the PSF magnitudes used for quasars). 

2.4. NVSS Catalog 

NVSS is a 1.4 GHz continuum sky survey with the Very 
Large Array (VLA) (Condon et al. 1998). It covers the sky 
region that lies north of 8 = —40°. The entire survey area is 
about 82% of the celestial sphere. The survey catalog includes 
about 2x10^ discrete objects with a full width half maximum 
of 45" and a nearly uniform sensitivity. The error in RA and 
dec varies from < 1" for sources with fluxes above 15 mJy to 
about 7" for sources at the survey flux limit of 2 mJy (Condon 
et al. 1998). The median redshift for the NVSS objects is 
around z = 1 . 

3. THE CROSS-CORRELATION AND SYSTEMATICS 

We have cross-correlated the SDSS QSO catalog with the 
WMAP sky maps to estimate the average amplitude of the SZ 
distortion associated with a known QSO. The estimated am- 
plitude of the quasar signal in a given WMAP band is simply 

AT = ^NiT}j ^N^j , (1) 

where A',' is the number of quasars in the ith WMAP pixel and 
Ti is the effective thermodynamic temperature of the pixel in 
the chosen frequency band. This expression just gives the av- 
erage temperature distortion associated with each individual 
QSO. The theoretical SZ signal from a single QSO is ex- 
pected to be about 1-10 juK for an arcminute beam depend- 
ing on the mass of the black hole (Scannapieco, Thacker, & 
Couchman 2008; Chatterjee et al. 2008). This when trans- 
lated to WMAP beams gives a signal of roughly 10-100 /Jy 
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Fig. 1 . — Sky coverage of SDSS (top row) and NRAO VLA Sky Survey (NVSS) (bottom row) in MoUweide projection. The right column shows the same 
maps after applying the radio (2 mjy) mask (SDSS only) and dust {E{B — V)> 0.05) mask (both). 



at WMAP frequencies (Refreiger, Spergel, & Herbig 2000). 
At WMAP frequencies, the SZ effect reduces the tempera- 
ture in a pixel, so the average temperature in Eq. 1 will be 
negative if dominated by the SZ effect. Galactic foregrounds 
and other systematic emission (described later) will positively 
bias the signal and lead to an underestimate of the SZ temper- 
ature decrement. For comparison, we also cross-correlate the 
SDSS LRG catalog with the WMAP sky maps, which serves 
as a rough measure of the galaxy cluster SZ signal to the ex- 
tent that the LRG sample predominantly traces clusters. 

Signals contributing to the sum in Eq. 1 are affected by a 
number of possible systematic effects, which we attempt to 
minimize. Sky regions with higher dust extinction will lead 
to selection biases in the QSO and LRG samples. Also, dust 
emission in these bands will contaminate the WMAP signal. 
We employ the dust prescription described in Ho et al. (2008). 
Using the Schlegel, Finkbeiner, & Davis (1998) extinction 
map, we construct E(B-V) masks excluding pixels with E(B- 
V) above some threshold value (see Table 3 for the dust mask 
thresholds). We apply the E(B-V) mask to both our SDSS and 
NVSS samples. Figure 1 shows the SDSS and NVSS sample 
spaces after applying the dust masks. 

Another source of contamination in the signal will come 
from radio emission from radio-loud quasars. The spectral 
index of these radio-loud sources is commonly around 0.7 
(Carlstrom, Holder & Reese 2002), so radio contamination 
at low frequencies is an important issue. We have matched 
the objects in the SDSS quasar catalog with the NVSS radio 
catalogs to find radio counterparts of the Sloan quasars. This 
was accomplished by (note that WMAP sky pixels are 7 ar- 
cminutes) choosing radio sources which are within an angular 
distance of 45 arcseconds, the beam size of the radio survey 
(although the pointing error is substantially less). We then 



mask WMAP pixels that contain selected radio sources with 
fluxes above 2 mJy at 1.5 GHz. The right column of Figure 1 
shows the Mollweide projections of our data sets after apply- 
ing the dust mask and the radio mask (for the SDSS samples 
only). 

We also want to minimize the noise arising from random 
associations of QSO positions with temperature distortions in 
the primary CMB anisotropics. Following the WMAP team's 
spatial Weiner filter analysis when obtaining point source am- 
plitudes from the temperature maps (Hinshaw et al. 2007), 
we use a similar high-pass spatial filter to cut off power at 
large angular scales where primary CMB fluctuations domi- 
nate the signal. Microwave temperature maps are convention- 
ally expressed as spherical harmonic coefficients, T{Q,(^) = 
£/„,a/„,y/m(9,c|)). A filtered map is reconstructed from the co- 
efficients flj^ = fiaim, with the Wiener filter (Tegmark & de 
Oliveira-Costa 1998; also see Chatterjee 2009 for detailed im- 
plementation of the filter) 

_ / A%Wi \ (y (2/ + l)wA ' 
^' [c,W^ + nJ[^ c,W^ + n, ) ' 

where c/ are the multipole moments for the primary CMB 
power spectrum, Wi are the moments of the angular WMAP 
beam profile, and «/ give the detector noise at each multipole 
moment. We computed the c/ values from the best-fit WMAP 
cosmological model using the CAMB code (Lewis, Challi- 
nor & Lasenby 2000) (although the filter is relatively insen- 
sitive to the precise set of multipole moments used), the de- 
tector noise n/ from the noise-per-differencing-assembly val- 
ues for WMAP (WMAP explanatory supplement: Limon et 
al. 2009), and the WMAP beam moments from Page et al. 
(2003). The maximum I values that we have used for con- 
structing the filters for different frequencies depend on the 
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TABLE 2 
Fits for the Spectrum 



Sample 


Model 


Fits for y 

10-7 


f 


BIC 


QSO 


SZ 


,V -1.611 .6 


2.58 


5.28 


QSO 


SZ+Dust 


>' = 7.0±3.4 


2.45 


4.65 


LRG 


SZ 


3? = 0.02 ±0.7 


2.75 


6.59 


LRG 


SZ+Dust 


>> = 5.3 ±2.45 


2.78 


4.98 




100 



Fig. 2. — The cross-correlation spectrum of the QSOs. The dashed line 
is the best-fit SZ spectrum, while the solid line is the best two-component 
fit using both SZ and dust. The best-fit values for the y parameter and the 
corresponding reduced values are listed in Table 2. 



angular resolutions of the corresponding bands: / = 1000 (Q- 
band), / = 1500 (V-band), and / = 2000 (W-band). The fil- 
ter suppresses power at scales which are dominated by the 
primary CMB temperature fluctuations or detector noise; the 
precise form of the filter in Eq. 2 gives a better signal-to-noise 
estimate of the amphtude of a point source superimposed on 
the primary CMB fluctuations and observed with the WMAP 
beam and detector noise. At WMAP angular resolutions, the 
SZ signal from an individual QSO or galaxy cluster can be 
treated as a point source. 

4. RESULTS 

Figure 2 shows the average radiation intensity difference 
per QSO in three frequency bands, obtained from Eq. 1. For 
clarity across bands, the thermodynamic effective tempera- 
ture in each band is converted to intensity. Figure 3 shows 
the same signal, except for the cross-correlation with LRG's 
instead of QSO's. The dashed lines shows the best-fit SZ dis- 
tortion with the usual spectrum 



A/ 



[xco\h{x/2)-A]y, 



(3) 



where x = hv/{KBTo). The soUd lines are a two-component 
fit incorporating both SZ and dust emission with a spectrum 
I{x) = Ax"*, corresponding to an effective dust temperature 
power law index of 2 (Gold et al. 2009). The best-fit y val- 
ues for all cases are shown in Table 2. The errors on the fit 
parameters are generated in the following way. We resimulate 
the data by drawing an error from a Gaussian distribution with 
root-mean-squared values given by the standard errors in the 
estimated signal for each band calculated from the variance in 



Ar and number of pixels used for averaging. We then perform 
the fits on the mock data and determine the standard deviation 
of the fit parameters among all mock datasets. The value of 
this standard deviation is provided as the estimate of error on 
the parameter. The values for the reduced and the Bayesian 
Information Criterion (BIC = X^ + ^log(n); k: degrees of 
freedom; n: number of data points; Liddle, Mukherjee, & 
Parkinson 2006) for the various cases are also listed in Table 
2. The information criterion tests whether the improvement in 
the fits is worth the extra parameters. The two-component fit 
gives a lower value; the improvement to the BIC indicates 
that the improvement from adding an additional parameter is 
enough to outweigh the extra model freedom from the addi- 
tional parameter. This conclusion should perhaps be taken 
with a grain of salt since only three data points with large er- 
rors are being fit. Note also that the two-parameter fit gives 
different values of y: the larger SZ component can better fit 
the suppressed point at 41 GHz, and the resulting lower SZ 
signal at 94 GHz can then be compensated by the increased 
dust contribution. The overall spectra for the two possibilities 
are significantly different, and measurements with smaller er- 
rors and data points at higher frequencies can distinguish be- 
tween the two cases. 

To investigate the significance of the observed anticorrela- 
tion, we randomize the positions of the quasars and recalcu- 
late the values of Al for a thousand realizations. We find that 
in 0.1% (Q band), 1.3% (V band), and 6% (W band) of all 
realizations, the correlation amplitude is more negative than 
in the real sample. The E(B-V) > 0.05 mask has been used in 
this case. We are thus confident that we are observing a real 
average CMB temperature decrement associated with QSO's. 
We change the thresholds of the dust masks in calculating the 
correlation for the randomized case. The results do not change 
with the change in the thresholds of the dust masks. 

To study the possible impact of foreground emission or ab- 
sorption, we investigate the change in the cross-correlation 
amplitude when we change thresholds used for each mask. 
Table 3 shows the effect on the cross-correlation amplitude 
of changing the threshold of our dust mask. The signal in 
the Q and V bands varies significantly with the thresholds of 
the dust mask; we note that the cross-correlation amplitude 
increases in these bands (that is, it becomes more negative) 
when we use more restrictive masks. The signal in the W 
band is always consistent with noise. The fact that the cross- 
correlation signal does not converge with increasingly strin- 
gent dust map masks indicates that the observed signal is due 
to a combination of SZ and dust components. We speculate 
that some of the actual SZ decrement is filled in by dust emis- 
sion in regions with higher dust. Also since dust wiU not be 
correlated with the quasar positions, the mean dust signal con- 
tributions to the QSO and LRG sample should be comparable. 
To check this we examined the trend of the signal with the 
threshold of the dust mask using the randomized sample for 
both LRGs and QSOs. The results are shown in Table 4.0. 
The values for the LRG and QSOs are rougly consistent and 
also we see the correlation value to be high positive (com- 
pared to the other bands) in W band. This likely means that 
the y-parameters derived here for the two-component fit in- 
cluding dust are closer to the real actual y-parameters. Our 
cross-correlation signals are unchanged when the radio mask 
threshold is changed from 2 mJy to 4 mJy, so the signal likely 
has no significant radio contamination. This is further born 
out because the signal is unchanged if we increase the cov- 
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TABLES 

Effect of Varying the Dusk Mask 



Band 


E(B-V) > 0.01 


E(B-V) > 0.05 


E(B-V) >0.10 


Q 


-739±468 Jy/sr 


-278 ±120 Jy/sr 


-229± 114 Jy/sr 


V 


- 1004 ± 414 Jy/sr 


-136± 105 Jy/sr 


-124 ±97 Jy/sr 


w 


-30 ±616 Jy/sr 


-122± 160 Jy/sr 


-46 ±147 Jy/sr 



TABLE 4 

Effect of Varying the Dusk Mask with Randomized QSO and LRG positions 



Band 


E(B-V) > 0.01 


E(B-V) > 0.05 


E(B-V) > 0.10 


Q (QSO) 


0.12 Jy/sr 


0.50 Jy/sr 


0.84 Jy/sr 


V (QSO) 


0.19 Jy/sr 


0.85 Jy/sr 


1.23 Jy/sr 


W (QSO) 


1.6 Jy/sr 


8.2 Jy/sr 


10.9 Jy/sr 


Q (LRG) 


0.01 Jy/sr 


0.10 Jy/sr 


0.01 Jy/sr 


V (LRG) 


0.0 Jy/sr 


0.12 Jy/sr 


0.81 Jy/sr 


W (LRG) 


5.8 Jy/sr 


10.12 Jy/sr 


12.1 Jy/sr 



erage of the mask to include all WMAP pixels within 20 ar- 
cminutes of any radio source. 

One of the major possible contributors to the cross- 
correlation signal is an SZ distortion from galaxy clusters as- 
sociated with the QSOs. A number of authors have previously 
reported evidence of SZ distortion in WMAP temperature 
maps (e.g; Bennett et al. 2003; Myers et al. 2004; Hemandez- 
Monteagudo et al. 2006; Spergel et al. 2007; Diego & Par- 
tridge 2010), most likely from unresolved galaxy clusters. To 
estimate the effect from galaxy clusters, we cross-correlate 
SDSS LRGs, which serve as a tracer of galaxy clusters, with 
the same filtered WMAP maps and the same masks used for 
the QSOs. Errors are also calculated in the same way; the 
signal has about the same statistical significance as the QSO 
signal. The cross-correlation spectrum is shown in Fig. 3; 
the amplitude and shape of the distortion is very similar to 
the QSO case shown in Fig. 2. The similarity in the cross- 
correlation spectrum suggests that both the signals might be 
coming from the same SZ source population. 

Chatterjee et al. (2008) shows that the SZ signal from QSOs 
will be dependent on redshift due to the time evolution of 
quasar feedback. We perform this test by binning the quasars 
in two redshift shces (z < 1.5 and z > 1-5). The results do not 
show any significant evolution with redshift. 

5. DISCUSSION AND PROSPECTS 

The SZ effect serves as a viable probe for detecting hot gas 
in AGN environments, though the signal is small. We have 
presented a first step toward detecting this signal by cross- 
correlating optical quasar catalogs with microwave maps. At 
frequencies lower than the null frequency, the SZ signal is a 
decrement which makes it unique among other potential sig- 
nals at those frequencies. However, oversubtraction of fore- 
grounds can also lead to decrements. Wright et al. (2009), 
ffinshaw et al. (2007), and Chen & Wright (2008) report de- 
tection of negative fluxes associated with point sources in 
the WMAP 3-year and 5-year point source catalogs. In the 
present work we find a deficit in microwave flux associated 
wifli bofli QSO's and LRG's from SDSS, which suggests SZ, 
but the fits to the observed spectrum are not sufficiently con- 
straining to establish this conclusion definitively. 

Any SZ signal is likely some combination of the signal 
associated with gas heated by QSO's and hot gas in galaxy 
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Fig. 3. — The cross-correlation spectrum of the LRGs. 



clusters. The observed correlation of flux decrement with 
the LRGs is similar to that for QSOs. This suggests that 
the sources of the signals might be similar. LRGs are trac- 
ers of galaxy clusters, and hence the SZ effect that we ob- 
serve could be originating from galaxy clusters associated 
with QSOs. Diego & Partridge (2010) have analyzed WMAP 
data by stacking regions around known X-ray clusters, find- 
ing an SZ signal from galaxy clusters in the range of 10 to 
20 /iK, corresponding to tens of mJy. The mean fluxes we 
get are of the order of 1 mJy, while the theoretical mean SZ 
signal from quasars is in the range 10 to 100 /iJy, depending 
on the mass of the black hole (Chatterjee & Kosowsky 2007; 
Scannapieco, Thacker, & Couchman 2008). Our signal is at 
least an order of magnitude higher than the signal that we ex- 
pect from quasars and an order of magnitude less than what is 
expected from massive galaxy clusters. 

The median redshift of the quasars in our sample is 1.3 — 
1.4; quasars at these redshifts tend to reside in halos with 
masses around a few lO'^M© (e.g.. Coil et al. 2007) and 
not in cluster mass haloes. However, since the cluster SZ 
signal is so much larger than the QSO SZ signal, our re- 
sults can be explained if only a small fraction of QSO's re- 
side in galaxy clusters. If 10% of QSOs reside in clusters 
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with 10^^ to 10'^M(:j halos, they would provide the observed 
cross-correlation signal. We made an estimate of the clus- 
ter signal from a typical halo mass for QSOs IO'^Mq and a 
typical cluster mass IO'^Mq by estimating the angular cross 
correlation function Tlgso.ciuster{Q), which will be equal to bias 
(QSOs) X bias (lO'^ solar mass halos) x mean dark matter an- 
gular correlation function at z « 1 . To get an estimate of the 
SZ contribution of the cluster signal we need to integrate the 
correlation function over the beam size which in Jy/Sr will 

be,(neg,„„,)"'/o''''™"'a'e27ien,,o,,,„,,,,(e)x(mean cluster SZ 
decrement) X (angular number density of clusters). Assum- 
ing a beam size of 20 arcminutes for WMAP, bgso = 2.31 
(Shen et al. 2009), bcu,ster = 3.3 (Estrada et al. 2009), a cluster 
density of 10 per square degrees (Carlstrom, Holder & Reese 
2002), the angular correlation function estimate at redshift 
1.0 (Abell et al. 2009), and a typical SZ signal from a clus- 
ter to be ^ 2QmJy we get an estimate of the cross-correlation 
to be ^ 65 Jy/Sr. This number is roughly comparable to our 
estimates of the cross-correlation. However this estimate is 
slightly lower than our detected signal which might imply 
that some of the power in the signal could be still related 
to backgrounds. An alternate possibility is that the SZ sig- 
nal from quasar feedback is higher (Natarajan & Sigurdsson 
1999) than those obtained from current simulations (Chatter- 
jee et al. 2008; Scannapieco, Thacker & Couchman 2008). 

Higher-resolution CMB maps can distinguish between 
QSO and galaxy cluster SZ signals by simply resolving the 
cluster signals (to all redshifts) and removing these pixels 
from any cross-correlation. Experiments which may improve 
on the current cross-correlation analysis include the Planck 
satellite, with angular resolution down to a few arcminutes in 
its higher frequency bands, and the ground-based arcminute 
resolution experiments like ACT, SPT, and APEX. The SZ ef- 
fect is subdominant at angular scales probed by WMAP (e.g., 
Komatsu & Seljak 2002; Bennett et al. 2003), which will not 
be the case at higher resolution. The mean QSO SZ am- 
plitude is predicted to be on the order of 1 /jK at 150 GHz 
for an arcminute size beam (Scannapieco, Thacker & Couch- 
man 2008), with substantial modeling uncertainty. Data from 
SDSS reveal around 30 photometrically detected quasars per 
square degree. 400 square degrees of sky coverage can pro- 
vide cross-correlation with around 10000 QSO's. If a ground- 
based experiment can attain a noise of 10 /jK per pixel over 
this sky area, this will give roughly a lOo detection of the 
QSO cross-correlation signal and a clear separation of the 
cluster and QSO SZ signals. Planck will have 7-arcminute 
resolution at 143 GHz; this implies a mean QSO SZ decre- 
ment of around 40 nK per pixel. The noise per pixel af- 
ter 14 months is projected to be around 6 /jK. With around 
5 X 10^ quasars over the full sky, the signal-to-noise of the 
cross-correlation will be 5a. Note that both of these estimates 
are for single frequency bands at the maximum decrement. 
Planck and ACT in particular are both mapping at frequencies 
above the null with higher angular resolution, giving stronger 
statistical constraints but more difficult systematic errors due 
to possible confusion between the SZ increment and fore- 
ground emission. Planck's many frequency bands may allow 
efficient separation of SZ from foreground through spectrum 
information. 

We have demonstrated a statistically significant correlation 
of a temperature decrement in WMAP sky pixels with the 
number density of both quasars and luminous red galaxies 
drawn from the Sloan Digital Sky Survey. The inferred QSO 



y-distortion for a two-component fit including both dust and 
SZ gives y = (7.0 ±3.4) x 10^^. The signal appears to be 
partly correlated with local dust emission. Higher signal- 
to-noise measurements are required to separate the dust and 
SZ contributions, while higher spatial resolution will sepa- 
rate QSO and galaxy cluster SZ signals. Both prospects are 
likely in the near term with the Planck satellite and dedicated 
ground-based SZ measurements. We anticipate that SZ mea- 
surements will soon make a valuable contribution to char- 
acterizing the energy feedback from quasars, to complement 
current data from the X-ray and radio bands. 
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